Spontaneous spin stripe dimerization in the doped t-J model. 
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To investigate spin dimerization in the t— J model we consider the extended t — J — 5 model with 
explicit spin dimerization introduced via parameter S. At zero doping the dimerized spin liquid is 
unstable at 5 < 5^ ~ 0.3. We demonstrate that the doping stabilizes the dimerized stripe phase: 
the Sc{x) decreases when doping x increases. At doping larger than critical, x > Xc, the dimerized 
phase is stable even without explicit dimerization, i.e. a.t S — 0. 
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It is widely believed that the 2D t — J model is rele- 
vant to the low energy physics of high-temperature su- 
perconductors. This is why investigation of this model 
is of great interest both for theory and experiment. In 
spite of great efforts during more than a decade there is 
no full understanding of the phase diagram of the t — J 
model, however some facts are well established. At zero 
doping the model is equivalent to the Heisenberg model 
on a square lattice which has long range Neel order 
Doping by holes destroys the order. A simplified picture 
of noninteracting holes leads to the Neel state instabil- 
ity with respect to spirals at arbitrary small but finite 
doping jsj. However more sophisticated numerical cal- 
culations which take into account renormalization of the 
hole Green's function under the doping indicate that the 
Neel order is stable below some critical hole concentra- 
tion Xc In the Neel phase [x < Xc), in all waves 
except s-wave, there is magnon mediated superconduct- 
ing pairing interaction between holes. Q. It is also clear 
that at very small hopping there is phase separation in 
the model because separation leads to reduction of the 
number of destroyed antiferromagnetic links |^ 

The purpose of the present work is to elucidate spin 
structure of the ground state in quantum disordered 
regime, i.e at a: > Xc- The most important hint comes 
from experiment: indications of stripes in the high- Tc 
materials 0. Another important hint is a remarkable 
stability of the spin dimerized phase in the frustrated 
Ji — J2 model. The idea of such state for this model was 
first formulated by Read and Sachdev and was then 
confirmed by further work |^,|l0|. The stability of such a 
configuration implies that the lattice symmetry is spon- 
taneously broken and the ground state is four-fold degen- 
erate. Such a route towards quantum disorder is known 
rigorously to take place in one dimension, where the Lieb- 
Schultz-Mattis (LSM) theorem guarantees that a gapped 
phase always breaks the translational symmetry. Some 
time ago Affleck suggested that the LSM theorem can be 
extended to higher dimensions, and the gapped states of 
quantum systems necessarily break the discrete symme- 
tries of the lattice |l|] . The example of the Ji — J2 model 
provides further support for this idea. 

There have been several attempts to consider the spin- 
dimerized phase in a doped Heisenberg antiferromag- 



net. For this purpose Affleck and Marston ana- 
lyzed Hubbard-Heisenberg model in the weak-coupling 
regime, Grilli, Castellani and G. Kotliar ||l^ considered 
SU{N), N — > 00, t — J model, and very recently Vo- 
jta and Sachdev (l4| considered Sp{2N), N ^ 00, t — J 
model with long range Coulomb interaction. These works 
indicated a stability of the spin-dimerized phase in some 
region of parameters, providing a very important guid- 
ing line. However relevance of these results to "physical 
regime" of the t — J model remained unclear. In the 
present work we demonstrate that the spin-dimer order 
is stable, the only small parameter used in the analysis 
is hole concentration with respect to the half filling. 
The Hamiltonian under consideration is 

H ^~tY^ cj^cja + (s,Sj - j . (1) 

is the creation operator of an electron with spin a 
{a =t: i) &t site i of the two-dimensional square lattice. 
The cj^ operators act in the Hilbert space with no double 
electron occupancy. The (ij) represents nearest neigh- 
bour sites. The spin operator is — 5 /3 '^L'''a/3Ci/3, 

and the number density operator is rii = J^a cLcio-. An- 
tiferromagnetic interactions Jy > are arranged in a 
stripe pattern shown in Fig. 1: solid links correspond 
to Jij = J± = J(l + S), and dashed links correspond 
to J,j = j = J(l - S). At half filling {{n,) = 1) the 
Hamiltonian ([|) has already been studied ^,|l^: for 
S > 6c ~ 0.303 the ground state is a quantum disor- 
dered state with gapped spectrum, and for S < Sc there 
is spontaneous Neel ordering with gapless spin waves. 

In order to study the stability of the dimer phase we 
first derive an effective Hamiltonian in terms of bosonic 
operators creating spin- wave triplets (magnons) tj^, a = 
X, y, z and fermionic operators creating holes a\, a =t, i 
from the spin singlets shown in Fig. 1. This Hamil- 
tonian consists of four parts: the spin- wave part Ht, the 
hole part Hi^ the spin- wave-hole interaction Hth, and the 
hole- hole interaction Hhh- Let us start from Ht- Similar 
effective theories have been derived in Refs. and we 
only present the result: Ht = H2 + H3 + H4 + Hu, where 

H2 = Ek,. {^k^L^ka + ^ {tljUo. + h-C.) }, H, - 
El+2=3R'(kl,k2)ea/37ikia*L/3*k37 + h.C, UlidHi = 
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We also introduce an infinite repulsion on each site, in 
order to enforce the kinematic constraint tj^tl^ — 0. 
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The following definitions are used in H2 '. -Bk ~ J (cos ky — 
O.Scosfca;) and = J± + -Bk- The matrix elements 
in the quartic and cubic interaction terms are: T(k) — 
j (0.25 cos /ci: + 0.5cosfcy) and R(p,q) = 0.25j(sin g^; — 
sin Pa;). Throughout the paper we work in the Brillouin 
zone of the dimerized lattice. 

At zero doping {{rii) = 1) iJj is an exact mapping of 
the original Hamiltonian ([^). To analyse this case it is 
enough to apply the technique [ p7pC| ] . The result for the 
normal spin-wave Green's function reads: 



GAr(k,tj) 
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{u + Ak{^uj)}{Lu - Akiuj)} + Bl 
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where ik(w) = ^k + (k) + E^^^(k,w) and Bk(w) 
-Bk + Sf (k). Normal S4' and anomalous self-energies 
are caused by the quartic interaction and the most 
important contribution S comes from the Brueckner 
diagrams as described in |17|. Strictly speaking there is 
also some contribution to the self-energy caused by the 
"triple" interaction -ffs. However this contribution is very 
small (see, e.g. Ref. [|lO|) and therefore we neglect it. 

Expansion of the self-energy in powers of w near 
Lo ^ gives quasiparticle residue and spin- wave spectrum 

^k ' = 1 - '^k = ^k\/[Ak(0)]2-Bk'. Expressions 

for effective Bogoliubov parameters Uk and are given 
in The spin- wave gap A = Wko, kg = (0,7r), ob- 

tained as a result of a selfconsistent solution of Dyson's 
equations is plotted in Fig. 2 (line at a; = 0). The crit- 
ical value of the explicit dimerization (point where the 
gap vanishes) 5c = 0.298 is in agreement with results 
of series expansions |^ and quantum Monte Carlo sim- 
ulations [p^ . The validity of the Brueckner approxima- 
tion is justified by the smallness of the gas parameter 
f^t = J2aiAa^ia) ■ At the Critical point Ut = 0.13. 

Consider now doping by holes. On the single dimer |s) 
the hole can exist in symmetric and antisymmetric states. 
Because of hopping between the dimers there is mixing 
between these states, but the mixing is very small (few 
per cent) and can be neglected. The symmetric state has 
substantially lower energy and therefore only this state 
is populated at doping. The corresponding hole creation 
operator is defined in the following way: ^/2al,\s) = 



(c| ^ + c[ cr)|0), where 1 and 2 numerate the dimer sites. 
Bare hole dispersion can be found by calculating hopping 
matrix elements. This gives Hh = X](cp + consi)a|,^apCT, 
where the const is chosen in such a way that 



Ep — t{coapy + 0.5 cos -t- 1.5) 
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vanishes at the minimum: = 0, po = (tt, tt). 

The spin-wave-hole interaction Hth can be easily cal- 
culated in the way similar to that for doped spin-ladder 
Jl8| , p^ . This interaction consists of two parts. The first 
one is interaction of a hole and a magnon positioned at 
different dimers. This is a relatively weak interaction 
which can be neglected |^^. The second part, which 
gives the main effect, comes from the constraint that a 
hole and a magnon can not coexist at the same dimer: 
Aa^la ^ 0. To deal with this constraint we introduce, 
similarly to (^, an infinite repulsion 
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The exact hole-magnon scattering amplitude caused by 
this interaction can be found via Bethe-Salpeter equation 
shown in Fig. 3a. It is similar to that for magnon-magnon 
scattering The result is 




(6) 



where E and k is total energy and total momentum of 
the incoming particles. 

Let us denote the hole concentration by x ~ n/N, 
where n is number of holes, and N is number of sites. 
Hence on-site electron occupation number is (nj) = 1 — x. 
Concentration of holes in terms of the dimerized lattice 
is two times larger n/(0.5A^) = 2x, and this is the gas 
parameter of the magnon-hole Brueckner approximation. 
According to (^) the holes are concentrated in the pocket 
in the vicinity of po = (tTjTt). Therefore the magnon 
normal self-energy described by the diagram Fig. 3b is 



2xr(cj,k + po) 
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It is instructive to consider first the case which al- 
lows an analytical solution: J± '3> j, V^nx <C 1. Bare 
magnon dispersion in this case is ujk ~ J± + j (cos ky — 
0.5 cos kx) and hence the integrals in (HJ^) can be calcu- 
lated analytically with logarithmic accuracy. This gives 
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where Su) = [uj ~ uj^ + j{uji^ ~ uju.^)/{t + j)] / {t + j), 
fi = max(|5a;|, v^TTx), and 9{5uj) is a step func- 
tion. The magnon Green's function is G'(k, uf) = 

^oj — ojk — S^^(k, w)^ . For illustration the spectral 

function ImG{uj) at k = ko, t/j = 3 and different x 
is plotted in Fig. 4. There are several conclusions from 
formula (||) and Fig. 4: 1) doping pushes the spin-wave 
spectrum up, 2) the effect is increasing with hopping t, 3) 
finite width appears, 4) there is only a logarithmic depen- 
dence on the infrared cutoff. Let us stress the importance 
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of the point (4). It means that the effect is practically 
independent of the long-range dynamics. Moreover, near 
the critical point (A = 0) the situation is even better: 
the spin- wave spectrum is linear and even the logarithmic 
divergence disappears. Thus in the 2D case there is sepa- 
ration of scales which justifies Brueckner approximation. 
If we tried to apply the described approach to the ID 
case (say a doped spin ladder) we would get into trouble: 
power infrared divergence appears in Brueckner diagram 
and hence there is no justification for gas approximation. 
Let us also comment on the point (3) (width). There is 
a "triple" contribution to the magnon self-energy, Fig. 
3c. This is a long-range contribution which is much less 
important than the Brueckner one, and this is why we 
neglect it (cf. with Ref. jl^). However this diagram 
influences the width of the magnon spectral function. 

In the general case there are two contributions to the 
Brueckner self-energy: S^', which is due to the magnon- 

(2) 

magnon constraint, and which is due to the magnon- 
hole constraint. To find the spin wave spectrum one 
has to solve selfconsistently Dyson's equation for Green's 
fimction (||), as it is described in Ref. Results for 

the spin- wave " gap" A as a function of explicit dimeriza- 
tion S for different hole concentrations x and t/ J = 3 are 
plotted in Fig. 2. Strictly speaking at x ^ the A is not 
a gap because of the large decay width. What we plot is 
the position of the centre of gravity of the magnon spec- 
tral function. However at A — + the width vanishes, and 
therefore the critical regime is uniquely defined. When A 
is not small there is an interesting question of calculation 
of the exact shape of magnon spectral function. However 
to resolve this problem one needs to include long-range 
dynamics (diagram Fig. 3c). 

It is clear from Fig. 2 that at t/J — 3 and x > Xc ^ 
0.090 the "gap" remains finite even a,t S ~ 0. This is 
regime of spontaneous dimcrization. For t/J ^ 2 the 
critical concentration is Xc — 0.106, and for t/J — 1, 
Xc = 0.132. Thus the doping stabilizes the dimerized 
phase. The larger the hopping t, the stronger the effect 
of stabilization. (The same follows from eq. (^).) This 
statement is true only if t/J ^ 10. At t/J ^ 10 there 
is a crossover to quasiparticles with higher spin (hole- 
magnon bound states) [|9| which indicate transition to 
the Nagaoka regime. The small parameter of the Brueck- 
ner approximation is concentration of holes in the dimer- 
ized lattice: 2x. Therefore at t/J — 3 one should expect 
~ 20% accuracy in calculation of Xc- Note that the value 
of Xc is close to that found in Q from the Neel state. 

An important effect related to the stabilisation is sup- 
pression of spin wave quantum fluctuations by doping. 
At X — and S — Sc — 0.298 the density of spin fluctua- 
tions is nt ~ 0.13. Increasing of the hole concentration to 
a; = 0.1 {a,t t/J = 3 and 6 = 6c) gives nt « 0.02. At the 
critical point t/J = 3, 6 = 0, x = Xc = 0.09 the density 
is just nt ~ 0.07 (see also |20t]). 



The phase diagram of the t — J — 5 model at zero tem- 
perature is presented in Fig. 5 Because of the mobile holes 
the dimerized spin liquid is a conducting state. Stability 
of this state is a very robust effect because it is due to 
the high energy correlations (typical energy scale ~ 2t). 
There are also low energy (long range) effects with typical 
energy scale ~ 2tx which can lead to hole-hole pairing, 
small amplitude density waves etc. We do not consider 
these effects in the present work because they are sec- 
ondary with respect to the main one: spin dimcrization. 
It is also worth noting that one can introduce additional 
parameters t' and t" (next and next-next neighbour hop- 
ping) which influence single hole dispersion 1^ , shift po- 
sition of the minimum and change shape of the "Fermi 
surface" : stabilization is not sensitive to these details. 

Strictly speaking one can not exclude a possibility 
that at decreasing S, at some point there is a first order 
phase transition to an absolutely different phase. How- 
ever such a transition at zero temperature always has 
a precursor: some specific bound state. The examples 
are: hole-magnon bound state for transition to Nagaoka 
regime in the present model at very large t, multi-magnon 
bound states for a frustrated spin-ladder as a precursor 
of spinon deconfinement pl| ], and multi-magnon bound 
states in 2D Ji — J2-niodel at the 1st order transition 
point J2/J1 ~ 0.6 ||l^. In the present model we do not 
see any hint of such precursors and therefore we conclude 
that the first order phase transition is highly unlikely. 

We have discussed in detail the transition to the Neel 
state at small hole concentration x. It is also clear that 
at large x there is a 2nd order transition to normal Fermi 
liquid shown schematically in Fig. 5 by the dashed line. 
One can say that it happens when hole concentration in 
the dimerized lattice is equal to unity, 2a:: ~ 1. However 
this is only a crude estimate. Unfortunately we can not 
describe this part of the diagram more precisely because 
our approach assumes that 2a; <C 1. 

It has been shown recently that Heisenberg antiferro- 
magnet on triangular lattice is very close to the insta- 
bility with respect to spontaneous dimcrization |2^ . We 
would like to note that it is highly likely that in this 
model doping also stabilises the dimerized phase. 

In conclusion, using the dilute gas approximation we 
have analyzed the phase diagram of the t — J — 6 model 
and the stability of spin-dimerized phase. At doping 
larger than the critical one, a; > a;c ~ 0.09, the dimerized 
phase is stable even without explicit dimcrization. 

I thank V. N. Kotov and M. Yu. Kuchiev for stimu- 
lating discussions. 
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FIG. 3. (a) Bethe-Salpeter equation for hole-magnon scat- 
tering vertex F. Solid line corresponds to the hole and 
dashed line to the magnon. (b) Brueckner contribution to the 
magnon self energy, (c) "triple" contribution to the magnon 
self-energy. 



FIG. 4. Magnon spectral density at k = ko = (0, n) in the 
limit J± S> j, t/j — 3, and different hole concentrations. 



FIG. 5. Phase diagram of the t — J ~ 5 model {t/J = 3) in 
the plane doping (a;) - exphcit dimerization (S). 



FIG. 1. Stripe spin dimerization on square lattice. Solid 
links correspond to J± = J{1 + S), and dashed links corre- 
spond to j — J{1 — S). 



FIG. 2. The magnon "gap" (centre of gravity of spectral 
function) as a function of exphcit dimerization S for t/J = 3 
and different values of hole concentration x, ijii) — 1 — x. 
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